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Cylindrical  Wave  Reciprocity  Parameter 

Robert  J.  Bobber  and  Gerald  9v.  Sabin 
U .  S.  Navy  Vndenoater  Sound  Reference  Laboratory,  Orlando,  Florida 
,  (Received  September  19,  1960) 

A  cylindrical  wave  reciprocity  parametei;  Jta-43/pf)-(4?X)l:^for  use  in  a  standard  reciprocity  calibration 
when  the  waves  between  the  source  and  receiving  transducers  are  cylindrical  is  derived  from  both  four- 
terminal  network  theory  and  wave  propagation  theory.  Experimental  data  from  calibration  measurements 


on  underwater  sound  line  transducers  are  presented 


INTRODUCTION 

The  cylindrical  v>'ave  reciprocity  parameter  is  the 
parameter  that  should  be  used  in  a  standard 
reciprocity  calibration^  when  the  shape  of  the  waves 
befeween  the  source  and  the  receiving  transducers  is 
cylindrical. 

The  reciprocity  parameter  used  in  conventional  reci¬ 
procity  calibrations  of  electroacoustic  transducers  as 
developed  by  MacLean®  is  a  spherical  wave  reciprocity 
parameter,  although  it  is  not  usually  identified  spe¬ 
cifically  as  such.  The  transmitting  response  of  the  re¬ 
ciprocal  transducer  is  defined  in  terms  of  the  sound 
pressure  produced  in  a  spherical  wave;  that  is,  there 
is  an  inherent  assumption  that  when  the  reciprocal 
transducer  is  transmitting,  the  hydrophone  or  micro¬ 
phone  being  calibrated  is  in  a  spherically  divergent 
sound  field.  This  requirement  of  spherical  waves  should 
not  be  confused  with  the  additional  standard  require¬ 
ment  that  a  hydrophone  be  calibrated  in  a  plane  pro¬ 
gressive  wave  sound  field.®  The  two  apparently  con¬ 
tradictory  requirements  are  not  incompatible.  In  prac¬ 
tice,  the  necessary  condition  is  achieved  by  using  a 
transmitter-to-receiver  distance  that  is  large  enough 
so  that  (a)  the  receiver  is  in  a  spherically  divergent 
wave  or  “far  field,”  and  (b)  the  radius  of  curvature  of 
the  spherical  wave  is  large  compared  to  the  receiver 
size  so  that  the  sector  of  the  wave  intercepted  by  the 
receiver  is  essentially  plane. 

plane  wave  reciprocity  parameter  has  been  de¬ 
veloped  by  Simmons  and  Urick.^  It  is  used  in  calibra¬ 
tions  where  all  the  transducers  are  plane  piston  radia¬ 
tors  and  the  transmitter-to-receiver  distance  is  very 
short  so  that  the  hydrophone  is  in  the  near  field  of  the 
transmitter.  Pulsed-sound  measurement  techniques  are 
required  to  avoid  standing  waves  between  the  parallel 
transmitter  and  hydrophone  plane  faces. 

Just  as  the  use  of  the  spherical  wmve  reciprocity 
parameter  /,  reejuires  that  the  transducers  be  effec¬ 
tively  points,  and  the  use  of  the  plane  w'ave  parameter 

'  American  Standard  Procedures  for  Calibration  of  Electroacouslic 
Transducers,  Particularly  Thosr^  for  Use  in  Water  (American 
Standards  Associations,  Inc.,  New  A’’ork,  1957),  Z24.24,  Sec. 
3.3. 1.2. 

“W.  R.  MacLean,  J.  Acoust.  Soc.  Am.  12,  140  (1940). 
^American  Standard  Acoustical  Terminology  (American  Stand¬ 
ards  .AssoriaLion,  Inc.,  New  York,  1960),  Sl.l,  Sec.  71.3. 

*  B.  D.  Simmons  and  R.  J.  Urick,  J.  Acoust « Soc.  Am.  21,  633 
(1949). 


prove  the  validity  of  the  parameter. 

- ^ 

Jj,  requires  that  the  transducers  be  effectively  parallel 
planes,  the  use  of  the  cylindrical  Vave  parameter  Jc 
requires  that  the  transducers  be  effectively  parallel 
lines.  The  parameters  Jp,  Jc,  and  J,  are  thus  used  for 
wave  propagation  in  one,  two,  and  three  dimensions, 
respectively. 

In  practice,  Jc  can  be  useful  for  the  calibration  of 
long-line  hydrophones  without  requiring  a  long  test 
distance  or  spherical  wave  correction  factors. 

The  parameter  Jc  is  defined  in  the  usual  way  as  the  • 
ratio  of  voltage  receiving  sensitivity  to  transmitting 
current  response.  The  receiving  sensitivity  will  be  the 
standard  free-field  voltage  sensitivity.®  Whether  the 
transducer  is  effectively  a  point  on  a  spherical  wave, 
a  straight  line  on  a  cylindrical  wave,  or  a  plane  on  a 
plane  wave,  the  ratio  of  output  voltage  to  input  free- 
field  sound  pressure  will  be  the  same  in  each  case.  The 
transmitting  response,  however,  depends  on  tlie  kind 
of  wave  being  propagated  and  where  the  sound  pressure 
is  measured.  As  used  in  the  definition  of  Jc,  the  trans¬ 
mitting  current  response  is  defined  only  for  a  trans¬ 
ducer  that  can  transmit  cylindrical  waves  and  only 
in  terms  of  the  sound  pressure  in  the  cylindrical  waves. 

As  will  be  shown  later,  an  acoustic  line  in  a  free  field 
will  transmit  cylindrical  waves  to  a  maximum  distance 
approximately  equal  to  L^/\,  where  L  is  the  length  of 
the  line  and  X  is  the  wavelength.  If  the  line  is  n  wave¬ 
lengths  long,  the  maximum  distance  is  n-\  or  nL. 

Two  different  and  independent  methods  have  been 
used  to  determine  the  cylindrical  wave  reciprocity  pa¬ 
rameter.  The  first  method  is  a  “network”  analysis  in 
which  the  familiar  technique  of  treating  an  electro¬ 
acoustic  system  like  a  four-terminal  electrical  network 
is  used.  The  second  method  is  a  “wave”  analysis  in- 
vohdng  Fresnel  integrals  and  the  known  spherical  wave 
parameter. 

NETWORK  ANALYSIS 

Figure  1  is  a  schematic  drawing  of  an  electroacoustic 
system  with  dimensions  or  boundaries  such  that  the 
transmitted  waves  are  cylindrical  at  a  radial  distance 
R.  The  system,  assumed  to  be  linear,  passive,  and  re¬ 
ciprocal,  is  described  by  two  linear  equations : 

(1) 

E=ZuU+Z,2l,  (2) 
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where  E  is  the  voltage  across,  and  I  the  current  through, 
the  line  transducer;  p  is  the  sound  pressure  measured 
at  a  radial  distance  R  or  anywhere  on  the  dashed  line ; 
U  is  the  volume  velocity  emanating  from  the  dashed 
line.  If  the  dashed  line  is  thought  of  as  a  second  line 
transducer,  U  is  the  volume  velocity  being  produced 
by  the  second  transducer,  or,  in  a  negative  sense,  it  is 
the  volume  velocity  being  absorbed  by  the  second 
transducer.  When  waves  merely  pass  the  dashed  line, 
U  is  zero.  The  coefficients  Zu  and  Z22  are  acijustic 
and  electrical  impedances,  respectively;  Zu  and  Z21 
are  electroacoustic  transfer  impedances  that  are  equal 
in  magnitude  if  the  system  is  reciprocal.  This  hypo¬ 
thetical  system  js  similar  to  that^sed  by  Wathen- 
Dunn®  to  derive  ^ 

First,  drive  the  electroacoustic  system  electrically 
with  a  current  Ii.  Let  pi  be^he  sound  pressure  in  the 
cylindrical  waves  at  a  radial  distance  R  from  the 
center  of  the  line  transducer,  or  anywhere  on  the 
dashed  line.  Since  the  waves  are  being  neither  produced 
nor  absorbed  along  the  dashed  line,  TJ  is  zero,  and 
Eq.  (1)  becomes 

pt—ZiJi. 

If  R  is  taken  as  a  standard  distance,  Zii=pijli  is  the 
transmitting  current  response  of  the.  line  transducer, 
or  5.  Therefore, 


t 


Fig.  1.  Cylindrical  wave 
electroacoustic  system. 
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M/S=U,/Pi. 


©Since  Jc  is  defined  as  M/S,- 


J.^M/S=U,/p,.  (7) 

Thus,  is  the  acoustic  transfer  admittance  between 
the  source  strength  U2  and  the  free-field  pressure  pi 
produced  at  distance  R. 

To  evaluate  Ui/pi,  the  pressure  in  a  cylindrical  wave 
must  be  analyzed.  It  is  given  by®'^ 


(3),  p=AHii^'^^{kR)  =  A\_J^{kr)-jN,{kRy], 


Next,  drive  the-  electroacoustic  system  acoustically 
with  a  volume  velocity  U 2  emanating  from  the  dashed 
line,  and  with  the  electrical  terminals  of  the  transducer 
open-circuited,  so  that  I  is  zero  and  the  voltage  is  £2. 
Then  Eq.  (2)  becomes 

Ei=ZnUi 

or 

Zii^Ei/Vi.  (4) 

Let  Pi  be  the  free-field  pressure  at  the  transducer 
location  when  the  system  is  driven  acoustically.  That 
is,  pi  is  the  sound  pressure  produced  at  radial  distance 
R  by  the  volume  velocity  emanating  from  the  dashed 
line.  If  the  numerator  and  denominator  of  Eq.  (4)  are 
multiplied  by  pi,  the  result  is 

Zii=  (£2/  pi)  {pi/  U  2). 

The  ratio  £2/ Pi  is  the  standard  free-field  voltage  sensi¬ 
tivity  or  M.  Therefore, 

Zii=  M  {pi/Ui),  (5) 

Since  the  system  is  reciprocal^^^nd  the  sign  or  phase  of 
the  parameters  are  unimportant,  O  ^ 

e 

mZli  —  Zii.  (6) 

Q 

Substituting  Eqs.  (3)  and  (5)  into  (6)  produces 

_  S==M{pi/Ui) 

‘  W.  Wathen-Dunn,  J.  Acoust.  Soc.  Am.  21,  542  (1949). 


where  k  is  the  wave  number  27r/X,-  R  is  the  radial  dis¬ 
tance,  4  is  a  constant,  £0®'  is  a  Hankel  function  of 
zero  order  and  second  kind,  /o  is  a  Bessel  function  of 
zero  order  and  first  kind,  and  £0  is  a  Neumann  function 
of  zero  order  and  first  kind.  As  kR  becomes  large,®'^ 

H„(^^ikR) - >{2/irkR)K  ■  (8) 

Figure  2  is  a  plot  of  both  and  the  approximation. 


Fig.  2.  Hankel  function  and  its  approximation  for 
_  large  values  of  the  argument. 


'  Philip  M.  Morse,  Vibration  and  Sound  (McGraw-Hill  Book 
Company,  Inc.,  New  York,  1948),  2nd  ed.,  p.  298. 

^  Leo  L.  Beranek,  Acoustic  Measurements  (John  Wiley  &  Sons, 
Inc.,  New  York,  1949),  pp.  57-62. 


O 


R.  J.  BOBBER  AND  G.  A.  SABIN 


448 


a 

© 


For  kR>0.35  the  approximation  error  is  less  than  10%. 
For  kR>3.0  the  error  is  less  than  1%.  Thus,  Eq.  (8) 
is  a  good  approximation  if  R  is  greater  than  about 
X/2,  and  .i 

p~A(2/TkR)K  (9) 

Morse®  evaluates  the  constant  A  as 

®  =  TT®  fpau,  (10) 

where  /  is  frequency,  p  is  the  density  of  the  medium,  a 
is  the  radius  of  the  cylindrical  radiator,  and  u  is  the 
linear  velocity  of  the  cylindrical  mdiator.  The  radius  a 
is  assumed  to  be  small  compared  to  a  wavelength. 

ISTo  assumption  has  been  made  as  to  how  the  electro- 
acoustic  system  was  dri^n  acoustically  or  as  to  the 
source  of  Ui.  It  is  assumed  now  thfl;  the  system  is 
driven  by  an  imaginary  line  transducer  located  on  the 
dashed  line  of  Fig.  1  with  an  infinitely  small  radius  a 
and  a  surface  linear  velocity  u.  This  assumption  in  no 
way  affects  the  development  of  Eq.  (7). 

Substituting  (10)  in  (9)  produces 

p=-irafpu(K/R)K  (11) 

The  volume  velocity  t/s  is  related  to  the  linear  velocity 
u  by 

U2=2iraLu,  (12) 

where  L  is  the  length  of  the  imaginary  line  transducer, 
and  also  the  length  of  tlie  real  line  transducer. 
Substituting  (11)  and  (12)  in  (7)  produces 

Jo=  U'i/pi=-  (2vaLu)/[jafpu(\/Ry\ 

,  ■={2L/pf){R/\y={2L/pc){Rk)K 


Fig.  4.  Logarithmic  plot  of  the  magnitude  of  Ei(W)  and  its 
^symptotic  e.xpressions  as  a  function  of  distance. 


the  receiver  are  effectively  point  sources,  and  the 
spherical  wave  reciprocity  parameter  can  be  used  for 
reciprocity  calibrations.  For  shorter  test  distances,  the 
appropriate  reciprocity  mrameter  is  derived  in  terms 
of  the  pressure  at  the  rrceiving  line  as  a  function  of 
separation  distance.  Each  element  %  the  line  source 
radiates  to  each  element  of  the  line  receiver.  The  ele¬ 
ment  of  pressure  on  dxi  om  the  receiving  line  that  is  due 
to  dx\  on  the  line  source  is  ^ 

.  dp=P^{\/r)e~^^^dx\,  .  (13) 

where  Po  is  the  reference  level  and  k  is  the  wave  number 
2x/X.  The  pressure  at  x^.  is  computed  by  integrating 
over  L,  the  length  of  the  line  source : 


where  c  is  the  speed  of  sound. 

WAVE  ANALYSIS 


P(X2)  =  RoJ’ 


The  second  approach  to  the  problem  of  determining 
the  cylindrical  wave  reciprocity  parameter  involves 
analysis  of  the  transmission  of  sound  from  a  line  source 
to  a  line  receiver,  as  in  Fig.  3.  The  two  line  transducers 
are  of  equal  length,  placed  parallel  to  each  other  and 
separated  by  a  distance  R.  When  R  is  sufficiently  large, 
the  pressure  at  the  receiver  is  inveBsely  proportional 
to  the  distance.  At  this  large  distance  the  source  and 


Fig.  3.  Line  source  and  a  line  receiver,  both  of 
length  L,  separated  by  distance  R. 


Similarly,  the  average  free-field  pressure  (^(P))av  ob¬ 
served  by  the  line.:-receiver  is  confuted  by  integrating 
p(x2)  over  the  length  of  the  receiving  line  and  dividing 
by  this  length,  so  that 

/Zf/2  ^Z(/2 

I  {\/r)e~^^^dxidx2.  (15) 

L/2'^-£/2 

With  the  u^al  Fresnel  approximations,  Eq.  (15)  leads 
to 

{p{R)\^^P^{\/L)E2{W),  (16) 

where  E2{W)  is  a  complex  second  Fresnel  integral 
C2{W)- jS2{W),^V=2L{2\R)-K  Details  of  this  inte¬ 
gration  are  shown  in  the  Appendix.  Figure  4  is  a  plot 
of  the  magnitude  of  £2  as  a  function  of  distance  R  in 
units  of  L-/\.  At  long  distances,  the  magnitude  of  £2 
drops  approximately  6  db  as  the  distance  is  doubled; 
this  is,  of  course,  the'spherical  wave  region.  At  short 
distances,  the  drop  is  about  3  db  as  the  distance  is 
doubled ;  this  is  the  cylindrical  wave  region.  The  transi- 
tiorf  from  the  cylindrical  wave  region  to  the  spherical 
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wave  region  occurs  at  R=L-/\.  This  is,  of  course,  a 
gradual  transition. 

The  asymptotic  values  of  Eq.  (16),  as  is  shown  in  the 
Appendix,  are 

{p{R))aw~Po(\/Ry  for  short  distances,  (17) 

~Po  (7^/72)  for  long  distances,  (18) 

so  that  the  expression  for  the  transmitting  current  re¬ 
sponse  S{R)  is  • 

@ 

S'(R)  =  {p(R))av/I=  (Pn/I)(X/Ry  for  short  distances 
S(R)—{P(i/I)(L/R)  for  long  distances. 

At  long  distances,  the  ratio  of  the  free-field  voltage 
sensitivity  to  the  transmitting  current  response  is  J,, 
the  spherical  wave  reciprocity  parameter.  Similarly,  at 
short  distances,  this  .ratio  is  Jc,  the  C3^ndrical  wave 
■  reciprocity  parameter.  With  the  earlier  condition  that 
the  free-held  voltage  sensitivity  is  the  same  for  spherical 
and  cylindrical  wave  conditions,  the  following  expres¬ 
sions' can  be  written: 

J,=M/S’XR) 

J,^M/S{R) 

7o=/,[5(i2)/5'(72)]=/,E(i?X)-5=  ilL/pc)(RX)i, 

which  is  the  same  as  the  result  that  was  obtained  by 
the  network  analysis  method.  . 

The  condition  Js=Jc  corresponds  to  the  transition 
distance  J?=Z,VX. 

MEASUREMENTS 

A  series  of  measurements  was  made  to  support  ex¬ 
perimentally  the  validity  of  the  cylindrical  wave  pa¬ 
rameter.  Three  underwater  sound  line  transducers,  each 
84  cm  long,  Were  calibrated  from  5  to  25  kc,  using  both 
a  near-field  cylindrical  wave  reciprocity  technique  and 
a  conventional  far-field  comparison  with  a  standard. 
Figure  5  shows  the  results  of  both  calibrations  of  one  of 
the  transducers.  The  agreement  between  these  two 
calibrations  is  good  with  maximum  differences  of  about 
1  db.  One  possible  source  of  error  for  the  near-fiield 
calibration  is  the  determination  of  the  effective  length 
of  the  acoustic  line.  The  poor  boundary  conditions  for 
the  long  test  distances  required  for  the  far-field  calibra¬ 
tion  are  also  a  possible  source  of  error,  particularly  at 
the  lower  frequencies.  Data  from  the  calibrations  of 


1  , 

8 -100 -I 


FAR'FtELD  CONVENTIONAL 
CALIBRATION 


1 


NEAR'FIELD  CYLINDRICAL 
WAVE  RECIPR^ITY 
CALIBRATION^ 


FREQUENCY  IN  KILOCYCLES  PER  SECONO 


Fig.  5.  Frcc-fickl  vnllagc  sensitivity  of  a  line  hvflrophone. 


Fig.  6.  Mea^red  relative  sensitivity  of  a  line  hvdrophone 
at  5  kc  as  a  function  ^  projector-to-hydrophone  measurement 
distance.  ®  * 

e 

the  other  two  transducers  show  essentially  the  same 
agreement.  ^ 

A  series  of  reciprocity  calibrations  as  a  function  of 
test  distance  is  shown  in  Fig.  6.  This  is  relative  receiving 
sensitivity  of  one  transducer,  at  one  frequency.  The 
same  measurement  data  were  used  for  both  curves  in 
the  figureT'One  curve  was  computed  by  using  Jc  in  the 
reciprocity  calibration  equations,  the  other  curve  was 
computed  by  using  /j.  At  short  distances  (less  than 
2  m  for  this  case)  cylindrical  J  gives  consistent  and 
accurate  results;  at  long  distances  (greater  than  4  m) 
spherical  J  gives  essentially  the  same  result.  Between 
2  and  4  m,  the  region  of  transition  from  cylindrical  to 
spherical  waves,  both  parameters  yield  inaccurate  re¬ 
sults.  The  maximum  discrepancy,  0.7  db,  appears  at 
the  distance  L^/\,  or  2.36  m  in  this  case.  The  calculated 
value  of  this  discrepancy  is  the  square  root  of  the  ratio 
of  the  magnitude  of  Eq.  (16)  to  Eq.  (17)  in  the  cylindri¬ 
cal  wave  region  or  to  Eq.  (18)  in  the  spherical  wave 
region.  The  discrepancy  appears  in  Fig.  4  as  one-half 
the  decibel  difference  between  t.t'e  curve  and  the 
straight  line  approximations.  The  value  at  t^  cross¬ 
over  distance  D/\,  which  should  be  the  maximum 
value,  is  calculated  as  —0.^  db.  This  value  agrees 
well  with  the  experimental  data  of  Fig.  6.  The  expected 
discrepancies  of  lesser  magnitude  predicted  by  Fig.  4 
for  the  cyljjidricaTwave  region  were  not  observed  in 
this  set  of  experin^^tal  data.  The  crossover  point  is 
.  proportional  to  frequency  so  that  at  25  kc  it  would  be 
at  almost  12  m  thus  illustrating  the  l(#ig  distances 
nee^d  for  conventional  spherical  wave  reciprocity  cali- 
brauon  ^  lin^ransducers.  •  * 

Measurements  on  other  line  transducers  have  yielded 
equally  satisfactory  results  for  reciprocity  calibrations 
at  short  distances  by  the  use  of  cylindrical  /.  In  par- 
^  ticular,  one  stave  (4§  wavelengths  long  and  J  wave¬ 
length  wide)  of  a  sonar  transducer  was  calibrated  by 
self-reciprocity  impedance  measurements*  using  short 
test  distances.  The  free-field  voltage  sensitivity  com¬ 
puted  by  Jc  showed  good  agreement  with  the  sensitivity 
determined  from  conventional  far-field  measurements. 

Although  the  three  transducer  reciprocity  measure¬ 
ments  were  made  with  three  transducers  of  the  same 

®  O.  A  Sabin,  J.  Acoust.  Soc.  Am.  28,  705  (1956), 
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length,  this  is  not  a  requirement  of  the  calibration. 
The  length  L  used  in  the  theory  applies  specifically  to 
the  length  of  the  reciprocal  transducer.  The  transducer 
used  only  as  a  sound  source  may,  in  any  type  of  reci¬ 
procity  calibration,  have  any  size  or  configuration 
provided  it  insonifies  the  other  two  transducers  with 
the.  same  effectively  plane  .sound  field.  The  length  of 
the  third  transducer  or  hydrophone  may  be  smaller 
than  L,  and,  within  the  assumed  condition  of  uniform 
cylindrical  waves,  be  only  a  point;  however,  the  cali¬ 
bration  of  a.^short  line  or  point  hydrophone  by  a 
cylindrical  wave  reciprocity  method  offers  no  advan¬ 
tage  over  a^^,,  conventional  spherical  waven,  reciprocity 
calibration.  More  important,  the  cylindrical  wav^rom 
the  reciprocal  transducer  would  not,  in  pracMce,  be 
perfectly  uniform..  The  sound  pressure  at  a  single  point 
^uld  probably  deviate  from  the  theoretical  uniform 
pTessure  much  more  than^ould  the  average  pressure 
along 'll  line.  Consequently,  the  calibration  of  a  point 
hydrophone  would  Be  subject  to  larger%rrors  than  that 
of  a  long  line  hydrophone. 

' '  CONCLUSION  C- 

It  is  interesting  to  compare  the  spherical,  cylindrical, 
and  plane  wave  reciprocity  parameters : 

J,=  (2/pc)(iiX) 

Jc={2/pc){RK)^L 

/p=(2/pc)(i?X)M. 

The  factor  (2?X)°  is,  of  course,  superfluous,  and  is  not 
normally  included  in  the  expression  for  Jp.  It  has  been 
added  here  to  illustrate  that  the  power  of  RK  in  each 
expression  is  the  same  as  the  power  of  R  in  the  spread¬ 
ing  loss  function  for  each  kind  of  wave.  That  is,  the 
sound  pressure  spreading  loss  is  proportional  to  R  for 
spherical  waves,  to  R  for  cylindrical  waves,-  and  to 
R°=l  for  plane  waves. 

The  effective  size  of  the  transducer  is  finite  for 
cylindrical  waves  and  plane  waves,  and  the  finite  di¬ 
mensions  appear  as  L,  the  length  of  the  line  transducer, 
and  A,  the  area  of  the  plane  transducer. 

All  the  parameters  contain  the  constant  2/ pc. 

Thus,  not  only  has  the  cylindrical  wave  reciprocity 
parameter  been  developed  by  two  independent  methods 
aiM  been  verified  by  measurements,  but  it  is  consistent 
wufi  the  spheri'cal  and  plane  wave  parameters  as  well. 


APPENDIX 

If  tlm  usuS  Fresn^lipproximations,  r=R  for  mag¬ 
nitude,  and  r=i?-|-(.Vi— .T2)-/22?  for  phase,  are  made  in 
Eq.  (15),  thefbsult  is 

(piR)\.=  iP,/LR)  ®  ^ 


exp[]— (.r  1 — .Vo)-/  2R2d.Vid.Xo. 

« 


(19) 


This  integral  has  the  form 


If 


and 


then 


I=A  f  f  f{v-^)drid^. 
J*  f(x)dx=glx) 

/g(xfdx=h(x), 

(ft 


(20) 


I==A\^h(a-{-b)-{-h(—a—b)—hia—b)  —  h{—d-{-b)'].  (21) 
In  particular,  if  a=b,  then 


I=Alh{2a)-\-h(-2a)-2h{0)'].  (22) 

In  the  case  being  considered,  as  is  shown  below  by 
series  representation,  h  is  an  even  function  so  that 
h{2a)  =  h{—2a)  and  ^(0)  =  0.  (Then 

•  I=2Ah{2a).  (23) 


Table  I.  Computation  of  the  function  E2{W)  —  C2iW)—jS20V). 


w 

CAW) 

■SsflE) 

IEAW)I 

0.0 

0.0000 

0.0000 

0.0000 

0.1 

0.0050 

0.0000 

0.0050 

0.2 

0.0200 

0.0002 

0.0200 

0.3 

0.0450 

0.0010 

0.0450 

0.4 

0.0798 

0.0034 

0.0799 

0.5 

0.1243 

0.0081 

0.1246 

0.6  . 

0.1781 

0.0168 

0.1789 

0.7 

0.2403 

0.0308 

0.2422 

0.8 

0.3095 

0.0517 

0.3138 

0.9 

0.3841 

0.0811 

0.3926 

1.0 

0.4616 

0.1200 

0.4769 

1.1 

0.5390 

0.1687 

0.5648 

1.2 

0.6132 

0.2269 

0.6538 

1.3 

0.6812 

0.2926 

0.7414 

1.4 

0.7404 

0.362$ 

0.8245 

2i 

0.7480 

0.3731 

^.8359 

1.5 

0.7898 

0.43*1 

0.9011 

1.6 

oMol 

0.5010 

0.9696 

1.7 

0.8640 

0.5606 

1.0300 

1.8 

0.8966  ^ 

0.6105 

1.0847 

1.9 

0.9326 

0.6516 

1.1377 

2.0 

0.9766 

0.6868 

1.1940 

i  2.1 

1.0298 

0.7223 

1.2577 

-  2.2 

1.0913 

0.7632 

1.3317 

2.3 

1.1551 

0.8138 

1.4130 

2.4  O 

1.2148 

0.8730 

1.4960 

2.5 

1.2653 

0.9356 

1.5736 

2.6 

1.3071 

0.9945 

1.6424 

tf.7 

1.3459 

1.0446 

1.7037 

2.8 

1.3882 

1.0862 

T.7626 

2.9 

1.4398 

1.1258 

1.8277 

3.0 

1.4988 

1.1706 

1.9018 
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Then 


J  jkx^/2R)dx 

=  j'  cos(kx^/2Ii)dx—jJ'  sm{ky?/2K)dx 

=  alC{x/cL)-iS{x/a)-],  (24) 

where 

a=(Xi?/2)l 

cos(Tf/2)dt—  j  I  sin(x/V2)ft^ 

0  *^o 

a  I  [C(x/a)  —  JS(x/a)Jdx=a^[C2(x/a)  —  jS2(x/a)J  (25) 

C,(u)-jSiiu)^  (  C{t)dl-j  f  S(t)dt. 

S?  *^0  *^0 


By  using  Eq.  (23)^ with  Eq.  (19),  the  equivalent  of 
Eq.  (16)  is  obtained: 

(;^(E)).v=Po(X/E)CC2(L/a)-i52(E/a)].  ’  (26) 

The  expressions 

C2{u)  =  uC(u)—  (l/ir)  sinJxM- 
Si{u)  =  uS(u)+  (l/7r)(cos§irM^— 1) 


given  by  Xlagnus  and  Oberhettinger®  can  be  used  for 
computations  involving  the  function  in  Eq.  (26).  These 
expressions  lead  to  the  values  given  in  Table  I,  which 
were  computed  for  values  of  C{u)  and  S{u)  in  the 
range  m=0  to  m=3.0. 

For  large  values  of  u,  C^iu)  becomes  approximately 
m/2  and  similarly  S2{u)  becomes  approximately  m/2. 
Thus,  when  u  is  large,  \E2{u)\~u/'<j2,  which  leads  to 
the  result  expressed  in  Eq.  (17)  for  short  distances. 
For  small  values  of  u,  the  series  must  be  written  for 
C2(m)  and  52(m),  as  follows: 


C2{u)  =  \u^—  (|t) 


Ik  w”  u' 

2 - bdir)^ - 


2-31-5 


2-5I-9 


(-l)n-l(|^)2(n-.)_ 

Cl  f  • 


2(2»-1)!(4m-3) 


S2(uJ=i‘ir- 


-(hr 


+(hr 


w 


rlP 


(28) 


2-21-3  ■■  2-41-7 

+  ...(_  l).-l(l^)2n-l. 


2-61-11 

m‘" 


2(2m)!(4m-1) 


These  series  result  fr*^m  integrating  twice  the  cosine 
and  sine  series,  respectively,  for  the  argument  •  Ittm®. 
When  M  is  sufficiently  small,  |E2(m)  1  =  \‘^,  which  leads 
to  the  result  in  Eq.  (18)  for  long  distances. 


®  W.  Magnus  and  F.  Oberhettinger,  Formulas  and  Theorems  for 
the  Functions  of  Mathematical  Physics,  translated  from  the 
German  by  John  Wermer  (Chelsea  Publishing  Company,  New 
York,  1954),  p.  96. 
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